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ABSTRACT 



Context. The discovery of true solar analogues is fundamental for a better understanding of the Sun and of the solar system. Despite a number 
of efforts, this search has brought only to limited results among field stars. The open cluster M67 offers a unique opportunity to search for solar 
analogues because its chemical composition and age are very similar to those of the Sun. 

Aims. We analyze FLAMES spectra of a large number of M67 main sequence stars to identify solar analogues in this cluster. 
Methods. We first determine cluster members which are likely not binaries, by combining proper motions and radial velocity measurements. 
We concentrate our analysis on the determination of stellar effective temperature, using analyses of line-depth ratios and Ha wings, making a 
direct comparison with the solar spectrum obtained with the same instrument. We also compute the lithium abundance for all the stars. 
Results. Ten stars have both the temperature derived by line-depth ratios and Ha wings within 100 K from the Sun. From these stars we derive, 
assuming a cluster reddening E(B - V) = 0.041, the solar colour (B — V) e = 0.649 ± 0.016 and a cluster distance modulus of 9.63. Five stars 
are most similar (within 60 K) to the Sun and candidates to be true solar twins. These stars have also a low Li content, comparable to the 
photospheric abundance of the Sun, likely indicating a similar mixing evolution. 

Conclusions. We find several candidates for the best solar analogues ever. These stars are amenable to further spectroscopic investigations and 
planet search. The solar colours are determined with rather high accuracy with an independent method, as well as the cluster distance modulus. 



Key words, stars: fundamental parameters - open clusters and associations: individual: M67 - stars: late-type 



1. Introduction 

The specificity of the Sun and of our solar system have been 
the subject of active investigation over the last 5 decades. How 
typical is the Sun for a star of its age, mass, and chemical com- 
position? How typical is that solar-type stars host planetary sys- 
tems? Are they similar at all to ours? 

The quest to find stellar analogues to the Sun has been 
going on for a long tim e (for an extensive review see, e.g., 
Cavrel de Strobelll 9961. and it stems from the poor knowledge 
we have of the Sun when seen 'as a star' and from how typ- 
ical the Sun is for a G2 type star, for its age, chemical com- 
position, populatio n. It is, however after the discovery of the 
first exo-planets (Mayor & Oueloz 1995 ) that this quest be- 
came even more compelling, because to find stars similar to 
our own would allow us to answer to fundamental questions 
related to the origin of the solar system, the frequency of plan- 
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* Based on observations collected at the ESO VLT, Paranal 
Observatory, Chile, program 278.D-5027(A). 



etary systems similar to ours, and ev entually the formation o f 
life in other exo-planetary systems (ICavrel de Strobellll996j) . 
The need to identify in the night sky solar proxies to be used 
for spectroscopic comparison is also diffuse, in particular for 
the analysis of small solar system bodies (Bohnhardt, private 
communication). 

Among th e mos t recent results in this research, 
" d2006h used high resolution, high signal- 



Melendez et al 



to-noise ratio Keck s pectra to show th at HD 98618 is a very 
close solar twin, and Ki ng et al .] d2005l) proposed HD 143436 
after analyzing 4 stars pre-selected from literature. These stars 
seem to compare well with the best known solar twin , HR 



6060, first analyzed by iPorto de Mello & da Silva j 

subsequently confirmed by ISoubiran & Triaudl ( 2004 
made a comparati ve study of several hundred s of ELODIE 
spectra. Finally, iMelendez & Ramfrezl (12007b have shown 
HIP 56948 to be the best solar twin known to date both in 
stellar parameters and in chemical composition, including a 
low lithium abundance. 



2 



L. Pasquini et al.: True solar analogues in M67 



The open cluster M67 is a perfect targe t to search for so 



lar analogues. Recent chemical analyses (ITautvaisiene et al 
2000; Randich et all 120061: IPace et all 120081) . 



10 



show that this 

cluster has a chemical composition (not only Fe, but also 
all the other elements) extremely similar to the solar one, 
as close as allowed by the high precision of the mea- 
surements. The analysis resulted in [Fe/Hl =- 0.03+0.03 for 
Tautvaisiene et alJd2000h . [Fe/H]=0.03+0.01 forlRandich et all 
(2006), and [Fe/H]=0.03+0.03 for lPace et all d2008h . 



There are other two additional characteristics which make 
M67 strategical. The first one is that all the determinations of 
age give for t his cluster an age encompassing that of the Sun 
(3.5-4.8 Gvr; lYadav et al.ll2008l) . while the age determination 
for field stars is always uncertain. The second characteristic is 
that M67 is among the ver y few clusters showing Li depleted G 
stars (Pasquin i et al.ll 1997p. This is an impo rtant point because, 
as pointed out bv lCavrel de Strobell (11996), even if many stars 
appear to have most characteristics similar to the Sun, their Li 
abundance is usually 10 times higher than in our star. Since Li 
is likely an indicator of the complex interaction taking place in 
the past between the stellar external layers and the hotter inte- 
rior, the choice of stars which also share the same Li abundance 
with the Sun is an additional property to pinpoint the true ana- 
logues. 

In our opinion, the search of analogues to the Sun and to 
the solar system can be well performed in open clusters (OCs), 
which show a homogeneous age and chemical composition, 
common birth and early dynamical environment. As a conse- 
quence, they provide an excellent laboratory for investigating 
the physics of solar stars and of planetary system evolution, 
besides being excellent probes of the structure and evolution of 
the Galactic disk. 

M67 is a rich cluster, therefore it provides us with the op- 
portunity to find many stars candidates sharing similar charac- 
teristics, and not only one. This is fundamental to obtain some 
meaningful statistics, and the cluster hosts many main sequence 
(MS) stars of mass around the solar mass, which form a con- 
tinuous distribution (Fig.Q]). 

Finding several solar analogues in M67 will also help in 
providing an independent estimate of the solar colors, a quan- 
tity which still suffers of some relevant uncertainty (see, e.g., 
Holmberg et al. 2006), as well as an independent estimate of 
the distance modulus of the cluster. 

The present paper is the culmination of a work, 
whic h involved the chemical determina tion of this clus- 
ter dRandich et all 120061: IPace et al.1 120081) . photometry and 



astrometry (lYadav et al.1 120081) to obtain membership, and 
FLAMES/GIRAFFE high resolution spectroscopy to clean this 
sample from binaries, and to look for the best solar ana- 



logues u^in^_fhehnej;depth ratios method dGrav & Johanson 



199 it iBiazzo et al.1 12007b 



(Cayrel de Strobel & Bentolila 



and the wings of the 
1989) to determine 



Ho- line 
accurate 

temperatures with respect to the Sun. In addition, the Li line 
is used to separate Li-rich from Li-poor stars. 



-5 13 

O 

> 14 
1 5 



1 7 




.•to. •.. 

m amm 



0.4 0.6 0.8 1.0 

B-V 



.2 



Fig. 1. Portion of the colour-magnitude diagram of M67. Our 
selected targets encomp ass the solar colou r, and are high prob- 
ability proper motion dYadav et al. 20081) and radial velocity 
single members. The red points refer to the stars observed with 
FLAMES/GIRAFFE in three nights. 



2. Observations and data reduction 

We obtained 2.5 hours in three observing nights in service 
mode with the DDT program 278.D-5027(A). Observations 
were carried out with the multi-object FLA MES/GIRAFFE 
spectrograph at the UT2/Kueyen ESO-VLT dPasquini et al 



|2002|) in MEDUSA modfl we were able to observe 90 targets 
(Table [ATTb . The employed setting was HR15N with central 



wavelength 665.0 nm, which with a wavelength range between 
644 and 682 nm covers simultaneously the Ha and the Li 1 res- 
onance doublet at 670.8 nm with a resolution of R~17 000. 
Three separate exposures were obtained to be able to identify 
short and intermediate period binaries by comparing the radial 
velocities at different epochs; the combined spectra have a typ- 
ical signal-to-noise (S /N) ratio of 80-1 10/pixel. 

We selected from the catalog of lYadav et al. (2008) the 
main sequence stars (13™0 $ V $ 15™0) with B-V close to that 
of the Sun («0. 60-0. 75) with the best combination of proper 
motions parameters, that is a membership probability supe- 
rior to 60%, and exclusion of candidates with a proper motion 
larger than 6 mas/yr with respect to the average cluster mem- 
bers. Full details about proper m otion errors and se lection cri- 
teria can be found in the original Yadav et al. (2008) work. 

The log of the observations is given in Table Q] The obser- 
vations were reduced using the ESO-GIRAFFE pipeline. 



1 This is the observing mode in FLAMES in which 132 fibers with 
a projected diameter on the sky of 1'.'2 feed the GIRAFFE spectro- 
graph. Some fibers are set on the target stars and others on the sky 
background. 
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Table 1. Log of the observations. 
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Date 


UT 


fexp 


DIMM seeing 






(d/m/y) 


(h:m:s) 


(s) 


(arc sec) 


132.875 


11.833 


06/02/2007 


06:24:29 


2200 


0.6 


132.875 


11.833 


1 1/02/2007 


04:01:08 


2200 


1.1 


132.875 


11.833 


23/02/2007 


01:34:55 


2100 


0.9 



Radial velocities were measured using the IRAI0 package 
FXCOR, which cross-correlates the observed spectrum with 
a template. As a template we used a solar spectrum acquired 
with FLAMES/GIRAFFE. Finally, the heliocentric correction 
was applied. The typical error for our single measurement is 
around 0.4 km s _1 . The three spectra/star were finally co-added 
to perform the spectroscopic determination of temperature and 
lithium abundance (see Sections [3.2l and [3~3T >. 

We note that the GIRAFFE solar spectrum^, taken with the 
same setup of our observations, is used through this work for 
spectroscopic comparison with the stars and the synthetic spec- 
tra. The solar spectrum has been obtained by averaging most of 
the GIRAFFE spectra (some show clear flat field problems and 
have not been used) and it has a nominal S /N ratio above 400. 



3. Data analysis and membership 

3.1. Radial velocity 

Out of the 90 stars observed, all selected on the proper mo- 
tion and membership criteria given above, we found that 59 
of them are probable single radial velocity (RV) members. We 
have retained all the stars which show RV variations smaller 
than «1 km s _1 in the three exposures acquired and which have 
a mean velocity within 2 sigma (»T.8 km s _I ) from the median 
cluster RV. In Fig. [2] the histogram of the radial velocity distri- 
bution of these stars is shown, together with a Gaussian fit with 
< y rad >=32.90km s" 1 and a cr = 0.73 km s" 1 . In Table[A2]ihe 
RV values are listed for the stars of the final sample, while in 
Table IA.3"1 the values of the single RV measurements are given 
for the stars we discarded. 

In Fig. [3] we show the enlarged portion of the colour- 
magnitude diagram CMD containing the original sample; in 
this Figure the discarded and the retained stars are indicated 
with different colours. Many of the discarded stars tend to oc- 
cupy the brighter side of the main sequence, where binaries are 
indeed expected to be present. On the other hand, our proce- 
dure still leaves several stars which are apparently above the 
photometric main sequence. This is because the radial velocity 
measurements are not of superb quality and because the time 
span by the observations is of only 18 days. Long period bi- 
naries will not be discovered by our three radial velocity ob- 
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Fig. 2. Histogram of the radial velocity distribution of the 59 
single members selected in M67 (continuous line). A Gaussian 
fit to the member stars distribution is also displayed (dashed 
line). 
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Fig. 3. Zoom of the M67 CMD, centered on the targets ob- 
served. In blue the 59 retained single member candidates are 
shown, in red the stars discarded. As expected by binary con- 
tamination, most of the discarded stars lay on the bright side of 
the MS. 



servations. We shall see as seven stars clearly stand up also in 
the Magnitude - Temperature diagram (see Fig. |5]l and they 
are best candidates for binaries of similar mass. We have kept 
them in the sample, and we anticipate that their presence does 
not influence our analysis or conclusions. 



2 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of the Universities 
for Research in Astronomy, inc. (AURA) under cooperative agreement 
with the National Science Foundation. 



3.2. Effective temperature 

Given that our targets are on the main sequence of a cluster 
of solar metallicity and age, the critical astrophysical param- 
http://www.eso.org/observing/dfo/quality/GIRAFFE/pipeline/solar.hBli4rs for the selection of the best solar analogues is the ef- 
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fective temperature. We have used two spectroscopic methods 
to compute the stellar effective temperature: the line-depth ra- 
tios and the Ha wings. To calibrate these methods we have 
used a grid of synthetic spectra, computed with SYNTHE 
from a grid of ID LTE model atmospheres computed with 
version 9 of the ATLAS code Jkuruczlll993a| bl) in its Linux 
version dSbordone et all 12004 ISbordonell2005l) . All the mod- 
els have be en computed with the "N EW" Opacity Distribution 
Functions (Castelli & Kurucz 2003) which are based on so- 
lar abundances from Grevesse & Sauvall d 1 998b with 1 kms 1 
micro-turbulence, a mixing-length parameter q-mlt of 1.25 and 
no overshooting. The grid of synthetic spectra covers the tem- 
perature range 5450-6300 K with [Fe/H]=0, logg=4.4377, 
£=1 km s and was degraded to the resolution of the 
FLAMES/GIRAFFE spectra. We stress that for both methods 
these models are used to quantify the difference between the 
stellar spectra and the solar spectrum. Zero point shifts are 
most likely present, due, for instance, to limitations in the at- 
mospheric models or to not perfect treatment of the Ha lines. 
While these inaccuracies will reflect in a wrong temperature 
for the Sun, the difference between the stars and the Sun will 
be much less affected. 



Table 2. List of the six line pairs used to derive the stellar tem- 
perature. 



AJA2 


LDR 


X1/X2 


(A) 




(eV) 


6469.210/6456.380 


Fei/Fen 


4.84/3.90 


6498.937/6516.050 


Fei/Fen 


0.96/2.89 


6608.024/6597.557 


Fei/Fei 


2.28/4.80 


6608.024/6627.540 


Fei/Fei 


2.28/4.55 


6646.932/6627.540 


Fei/Fei 


2.61/4.55 


6646.932/6653.850 


Fei/Fei 


2.61/4.15 



3.2.1. LDR method 

It has been demonstrated that for stars with B - V — 0.4 -1.5 
line-depth ratios (LDRs) are a powerful temperature indica- 
tor, ^ip_ablejoj£solyejemp 

K dGrav & Johanson|[l99l ICatalano et alJl2002t iBiazzo et af 



20071) . Since our stars are within this B - V range, we have 



applied the LDR method to the members previously selected 
by radial velocity measurements (see Table \A.2i . To convert 
the line-depth ratios of our stars into effective temperature we 
need to calibrate a temperature scale for the measured LDRs. 
To this purpose we have considered an initial sample of about 
100 lines of iron group elements (which are usually temper- 
ature sensitive) present in the spectral range covered by our 
observations, from which we selected lines with the following 
characteristics: weak (to avoid saturation effects), sensitive to 
temperature variations, and at the same time well measurable in 
our spectra. The final selection contains six line pairs suitable 
to apply the LDR method; we measured them in the synthetic 
spectra and then derived an LDR-T e ff calibration for each pair. 
In Table [2] we list for these six line pairs the wavelength, the 
element and the excitation potential, as taken from the NIST0 
Atomic Spectra Database Lines, while in Fig.|4]we show an ex- 
ample of one LDR-r e jf calibration. The methods for the mea- 
surement of the line depth and t he re lated uncertainties a re de- 
scribed in ICatalano et al.1 ( 120021) and lBiazzo et alj 020071) . In a 
summary, the lowest seven points in the core of each measured 
line were fitted with a cubic spline and the minimum of this 
cubic polynomial was taken as the line depth. Given the lim- 
ited S IN ratio and resolution, the errors in each line depth are 
dominated by the uncertainty of the signal in the continuum. 
We have then measured the LDRs and derived for each line 
pair the temperature; in Table IA.2I the averaged values for all 



National Institute of Standards and Technology. 



Fig. 4. Example of LDR-r e ff calibration obtained with syn- 
thetic spectra. The line-depth ratio is between /16498.937 A and 
,16516.05 A. 



line pairs are given, where the associated uncertainty reflects 
the scatter obtained by the six measurements. 

With this method the effective temperature of the observed 
GIRAFFE solar spectrum results in 5792+27 K, i.e. 15 K 
higher than the synthetic one (5777 K is the theoretical effectiv e 
temperature of the solar atmosphere: IWilson & Hudsonlll991 ). 

We computed the temperature difference Ar LDR between 
the FLAMES/GIRAFFE targets and the Sun (as obtained from 
the six line-depth ratios and the summed spectra of the targets) 
as a f unction of th e de-reddened B — V colour (E(B - V) = 
0.041 : lTavlorl2007l) . The Ar LDR -(B- V) relationship for our 
targets is well described by a linear fit, which gives A!T LDR = 
(-3662.65 + 351.22) x (B - V)o + (2410.59 + 216.52) and an 
rms of 100 K. 

In Fig.[5]the temperature-magnitude diagram is shown. The 
two colours of the symbols are referred to stars with lower 
and higher presence of lithium. Seven stars clearly stand out 
of the main sequence, suggesting a parallel binary sequence. 
They most likely are long period binaries with components of 
similar mass not detected as RV variable by our observations, 
because of the limited time base of our observations. 
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log(Nu) < 2 
log(N Li ) > 2 
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Fig. 5. V vs. r e ff of the probable candidates. The stars have dif- 
ferent colours according to the LTE lithium abundance. Seven 
stars depart from the main sequence; they are most likely long 
period binaries which escaped the detection in our observa- 
tions. 



3.2.2. Ha line 

The wings of 
to temperature, 



the Ha 
like all 



line profile 
the B aimer 



are very sensitive 
lines, and depend 



only slightly on metallicity and gravity ( Cavrel et al.l 11985 



Fuhrmann et al.l 1 1 993b iBarklem et alj|2002l) . In particular, the 
spectral region in the range between 3 and 5 A from the 
Ha line center is a good effective temperature diagnos- 



tic dCayrel de Strobel & Bentolilal[l989l) . With respect to the 
higher members of the series, the Ha line has considerable less 
blending in the wings, making easier the placement of the con- 
tinuum. A further advantage, over other members of the series, 
is that it is rather insensitive to convection and in particular to 
the adopt ed mixing length pa rameter aMLT in ID model atmo- 
spheres (IFuhrmann et al.lll993l) . Thus, we selected this region 
as temperature indicator and for each star we have compared 
the Ha line profile outside the core in our real spectrum against 
the synthetic profile. The hydrogen line profiles are computed 
in SYNTHE by routine HPR 0F4, which, for resonan ce broad- 
ening uses essentially the i & Grieml (119651 119661) theory, 
and for Stark broadening it calls routine S0FBET, writt en by 
Deane Peterson, which uses essentially the theory of iGriem 
dl960l) . with modif ied parameters, so as to provide a good ap- 
proximation to the IVidal et alJ (11973b profiles (F. Castelli, pri- 
vate communication). For further d etails on the computations 
of h ydrogen lines in SYNTHE , see ICastelli & Kuruczl (l200lh 
and ICowlev & Castellil (12002b . The dominant broadening for 
Ha is resonance broadening, while Stark broadening becomes 
relevant for higher members of the series. In order to minimize 
the subjectiveness of the measurement, we have quantified the 
comparison between the synthetic profile and the observed one 
minimizing the rms of the subtraction. Continuum normaliza- 
tion is not easy for such a broad line, however, the fact of 



using a fiber instrument with a large coverage, minimize the 
subjectiveness of the process and makes it quite reproducible. 
Given the limited S /N ratio of the observations, it is however 
very difficult to provide a realistic estimate of the involved un- 
certainties. The systematic errors in the effective temperatur e 
obtained from the Ha wings is given by iGratton et al.l (1200 lh . 



and t he errors associated to the method have been discussed by, 
e.g., iBonifacio et all d2007h . where the dominant source of er- 
ror for echelle spectra has been identified in residuals in the 
correction of the blaze function. The GIRAFFE spectra are 
fiber-fed and the flat-field is obtained through the same opti- 
cal path as the stellar spectra, thus flat-fielding allows a better 
removal of the blaze function than it is possible for slit spectra. 
We have estimated for our stars an average error of ±100 K. 
With this method the effective temperature of the Sun results 
in 5717+100 K, i.e. 60 K lower than the solar real value. We 
note that the absolute temperature determined in this way de- 
pends critically on a number of assumptions in the model, and 
on the adopted broadening theory for Ha, and these produce a 
zero point shift of the Sun. The relative measurements, which 
are made with respect to the observed solar spectrum, are in- 
stead rather insensitive to all the assumptions used to build the 
synthetic profile. 

A linear fit well describes the relationship between the 
temperature difference AT Ha of the FLAMES/GIRAFFE tar- 
gets and the Sun as obtained from the Ha wings, and the de- 
reddened B-V colour: AT Ha = (-3811.8±283.39)x(B- V) + 
(2488 + 174.71) with an rms of 81 K. 

The calibrations Ar e jf vs. (B - V)o obtained with the two 
methods agree quite well, as shown in Fig.[7j they have slightly 
different slopes which produce a maximum difference at the red 
edge ((B - V)q = 0.7) of the sample of 40 K (Ha temperatures 
are cooler). These relationships can be used to calibrate stars 
with metallicity close to solar. Our LDR T e g vs. ( B - V)p rela 



tionsh ip has almost exactly the same slope of the lAlonso et al 
(119961) relations hip, but it is h o tter th an this by ~60 degrees. As 



a reference, the lAlonso et al. I dl996l) 's scale produce an effec- 
tive temperature for the Sun of 5730 K for a B - V = 0.63. 

3.3. Lithium 

Lithium is an important element because it is easily destroyed 
in stellar interiors, and its abundance indicates the amount of 
internal mixing in th e stars. Lithium in P op I old solar stars 
varies by a factor 10 dPasquini et alJI 19941) and M67 is one of 
the few clusters w hich clearly shows th is spread among other- 
wise similar stars dPasquini et al ,1 ll997l) . 

Equivalent widths (EWs) of the lithium line at A — 
670.7876 nm were computed using the IRAF task SPLOT; 
from measured EWs we de rived Li abunda n ces u sing the 
curves of growth (COGs) of lSoderblom etalj dl993l) . At the 
GIRAFFE resolution the Li i lines are blended with the Fe i A 
670.744 nm line, whose contribution to the lithium blend was 
subtracted using the empirical correction of the same authors. 
Lithium abundances were then corrected for the NLTE effects 



using the prescriptions of lCarlsson et alJ (1994). 
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Fig. [6] shows the lithium abundance for the summed spec- 
tra of the 59 targets as a function of the effective temperature, 
as derived by the LDR method. The filled symbols refer to the 
LTE abundance, while the empty ones represent the non-LTE 
abundance. The difference between LTE and non-LTE values 
is minor. The blue points are the solar twins (see Section |4). 
The position of the Sun is shown at log ,/V(Li) LTE = 0.84 and 
r e L ft DR = 5792 ± 27 K ( see SectionEHl The lithium abundance 
is listed in Table IA.2I Several solar twin candidates have Li 
abundances which are comparable with the Sun, whose value is 
0.84 as measured by us on the GIRAFFE spectru m (see above 
and 1 .0 as measured in high resolution solar atlas (IMuller et 
19751) . In most investigations the error associated to the ef- 
fective temperature is usually the dominant one (100 K corre- 
spond to about 0. 1 dex in log N(Li)), but in this case, since we 
have good T e ff determination, but limited resolution and S /N 
ratio, the uncertainty in the abundance associated to the equiv- 
alent width measurements is not negligible. The expected un- 
certainty in th e measured lith ium equivalent widths has been 
estimated from CayreJ(11988)'s formula: 



e), 
al. 



1.6 



°EW Li = 77717 ^FWHM x 8x 
S IN 

where S /N is the signal-to-noise ratio per pixel, FWHM is the 
full width of the line at half maximum, and 6x the pixel size. 
The predicted accuracy, o"ew u , is 3.0 mA for a typical S /N 
ratio of 80 and of 1.6 mA for a S /N ratio of 150. However, 
it should be noted that this formula neglects the uncertainty 
on the continuum placement. We estimate that, using homoge- 
neous procedures for the determination of the continuum and 
the line widths, the statistical error for the weak lithium line is 
of the order of 2-3 mA, depending on the S /N ratio of the co- 
added spectrum. This will correspond to an asymmetric error 
A log N(Li) ~^_qI dex for a star with a line as weak as the Sun, 
and of +0.04 dex for a star with a log Af(Li)=2.2. Since the line 
is weak, the error (in percentage) is inversely proportional to 
the line strength. Given the errors, all the stars with upper lim- 
its in our sample may have a Li comparable to the solar one. 

After the early works on M67, several Li surveys have been 
carried out of additional clu sters well sampling the age metal- 
licity space dRandich|[ 2008 and references therein): out of nine 
clusters older than the Hyades with available Li measurements, 
only two, besides M67, show a significant dispersion. The latter 
seems to be an exception, rather than a rule and its occurrence 
does not depend on age, nor on metallicity, nor on global clus- 
ter parameters. 

In this context, the novel result of our analysis and, in par- 
ticular, of the careful selection and cleaning of the sample as 
well as of the precise effective temperature determination, is 
that the large spread is clearly present only for stars cooler than 
~ 6000 K. Stars warmer than 6200 K seem to show a decay, 
probably indicating the red side of the "Li-gap", while stars in 
the 6000< T e s < 6200 K do not show any major scatter. 

It is now well ascertained on empirical grounds that, in or- 
der to explain the MS Li depletion in solar-type stars, an extra 
or non-standard mixing mechanism must be at work. No con- 
sensus so far has been found on the nature of this mechanism; 
it is nevertheless clear that, whatever this process is, it must 
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Fig. 6. Lithium abundance versus T^ R for the most probable 
single stars observed in M67. Filled circles: LTE abundances. 
Open circles: NLTE abundances. The blue points indicate the 
best solar analogues candidates. The position of the Sun is also 
shown at logJV(Li) LTE = 0.842 and T^ R = 5792 + 27 K. 



be driven by an additional stellar parameter besides mass and 
chemical composition. The presence of the Li spread indeed 
indicates that this parameter must vary from star to star, and 
that, depending on it, some stars (including the Sun) undergo 
a much more efficient mixing than others, while the absence 
of a dispersion for stars warmer than 6000 K suggests that this 
parameter is more uniform among F-type stars. 

Recent modeling have had some success in reproducing the 
solar Li abundance by using fairly complex models wh ich in- 
clude internal gravity waves dCharbonnel & Talonl 120051) . how- 
ever, those models are not able to reproduce the observed evo- 
lution of Li with age, a nd, in particula r, the "plateau" in Li 
abundances at old ages (Randich 2008). We are not aware of 
similar models for a grid of masses, to be compared to our ob- 
servations; since the number of possible parameters which in- 
fluence the Li evolution is very large (depth of convective zone, 
initial rotation, magnetic field, mass losses and torques just to 
mention a few), we cannot really predict at present why the 
extra mixing takes places at a given T e ff in M67 stars. 



4. Solar analogues 

With the aim to find the best solar analogues, we have com- 
pared Ar LDR to AT Ha (Fig.0. 

There are in our sample 10 stars for which both the 7^jj? R 
and J e H ff a are within 100 K of the solar values (285, 637, 1 101, 
1194, e i303, 1304, 1315, 1392, 1787, 2018). We will use these 
stars to find the best solar analogues and our best evaluation 
for the solar colours and the cluster distance. These stars are 
indicated in bold face in Table IA.21 

The average difference between these 10 stars and the solar 
is of -13 K (with a sigma of 60 K), while the average 
is =-9 K, with a sigma of 58 K. The average 
characteristics of these 10 stars should therefore well represent 
the solar values. 



T LDR 
1 eff 

1 eff 1 eff 
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Fig. 7. Ar LDR as a function of AT Ha for the 59 probable single 
member stars. 



The average B - V of the ten analogues is < B — V >= 0.692 
(cr = 0.020), their average magnitude is < V >= 14.583 mag 
(cr = 0.147), and the < V - I >= 0.754 (cr = 0.025). The 
B—V spread is larger than the formal errors in the photometry, 
indicating a possible real spread in the stellar characteristics. 
This is not surprising because, formally, these stars may span a 
range of up to 200 K in temperature. 

If we take our results of Table lA.ll two stars (637 and 1787) 
have both T e s determinations within 50 K from the solar val- 
ues, and three additional ones (285, 1101, 1194) within 60 K; 
these 5 stars are overall the closest to the Sun, with nominal 
effective temperatures derived with both methods differing less 
than 60 K from the solar one. Their average magnitude (14.557 
mag) is very similar to what found for the full subsample of 
10, as well as their average B - V colour (0.688) just 0.007 
magnitude bluer than the whole subsample. 

All the data in our possession indicate that some of these 
stars have a metallicity very close (within 0.03 dex; note that 
also their Li abundance is comparable; see previous Section) to 
the Sun, they have a very similar temperature (within 50 K), 
as well as a comparable age to the Sun, and they are true main 
sequence stars. To the best of our knowledge they are the best 
candidates in M67 to be the closest analogues to our star. 

In Figure[8]we compare the GIRAFFE spectrum of the Sun 
with the sum of the spectra of the 10 best stars analogues and of 
the 5 best analogues in a portion of the spectra which includes 
Ha and in another including the Li lines. The extremely small 
difference between the solar spectrum and these co-added spec- 
tra confirm quantitatively the very close resemblance of these 
stars to the Sun. 

At the request of the referee we performed also a direct 
comparison between the solar spectrum and the spectrum of 
our solar analogs. We used a^ 2 minimization using a Doppler 
shift and a re-adjustment of the continuum of the stars of M67 
as free parameters, in order to match the observed spectra to 
the observed GIRAFFE solar spectrum. The reduced^ 2 of the 
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Fig. 8. Comparison between the GIRAFFE solar spectrum and 
the spectra of the 10 and 5 best solar analogues. The Ha and 
lithium regions are shown in the upper and bottom panels, re- 
spectively. 



fit, or the associated probability, then provides a mean to rank 
the M67 stars. We restricted the comparison to a range of 10 
A centered on Ha. While in the fitting of the synthetic spectra 
the core of the line was excluded from the fitting range, it was 
here included. The LTE synthetic spectra fail to reproduce the 
core of Ha due to the presence of a chromosphere (absent in 
the model atmospheres employed by us) and to NLTE effects. 
Instead, the sought-for solar analogs must behave exactly like 
the Sun, including in the core of Ha. With this method the three 
M67 stars which are most similar to the Sun are the stars 1 194, 
1101 and 637. The result is thus very similar to what obtained 
by comparing the observed spectra to synthetic spectra, con- 
firming that the stars we selected are very similar to the Sun. 
We prefer the method based on synthetic spectra, since the di- 
rect comparison to the solar spectrum is affected by the noise 
present in the latter. 
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5. Solar Colour 

We would like finally to use the observed B—V colour of the so- 
lar analogues to derive in an independent way the B—V colour 
of the Sun, and this requires to evaluate the cluster reddening. 

The reddening towards M67 has been evaluated by many 
auth ors in t he las t 50 years, and a thorough discussion is given 
by iTavlorl d2007l) . M67 reddening is evaluated by this author 
in E(B - V) = 0. 041 ± 0.004, wh ich is accidentally the same 



value obtained by lAn et al.l ([2007) as the average point of the 
traditionally accepted range for the cluster. We will therefore 
adopt this value. This implies that the de-reddened colour for 
the average of our 10 solar analogues is (B - V) = 0.651. 

The value is in excellent agreement with what found by 
inverting the fit of all the stars using the r LDR , which would 
have predicted (B - V) = 0.659 and the value obtained by 
inverting the fit of r Htt , which would give (B - V) = 0.651. 

It is not simple to evaluate a realistic error estimate for this 
colour. This should include: the spread (0.020 mag) around our 
determination, the uncertainty in the cluster reddening, plus 
other systematics originating from stellar evolution and pho- 
tometry. We evaluate the evolutionary effects by investigating 
the expected variations of the solar color with age and metallic - 
ity by using evolutionary models. Photometric uncertainties are 
estimated by comparing Yadav et al. (2008) photometry with 
what obtained for this cluster b y other gr oups. 

We use the tracks from Girar di et al. (2000) for analyzing 
differential evolutionary effects. Because our stars have a sim- 
ilar effective temperature to the Sun, age has no influence: for 
stars younger than 5 Gyr of solar T e g , the visual absolute mag- 
nitudes and colours do not change in any appreciable way. If 
M67 were younger than the Sun, the only effect would be that 
the masses of our stars were higher, by about 1 %, than the solar 
one, but no difference is predicted in magnitude or colours. The 
other source of systematic uncertainty is the possibility that 
metallicity is not exactly solar. In this case, for a fixed effec- 
tive temperature, we do expect that a star more metal rich by 
0.05 dex would be slightly brighter (0.08 magnitude in V) and 
slightly redder (0.01 mag) than the Sun. 



Our phot ometry is taken by Yadav et al.l (120081) . which was 
calibrated on lSandquist et al. ( 2004 ). In order to check for pho- 
tometric systematic errors, we have compared our valu es for 



the 10 best analogues with Mo ntgomery et alj (119931) . find- 
ing that using their photometry an average (B - V)o=0.650 
would have been foun d, i.e. only 1 mmag bluer than our value. 
Sandquist et alj d2004l) made on the other h and a general com- 



the M67 stars to compensate for their higher metallicity. We 
conclude that a solar (B - V) = 0.649 + 0.016 is our present 
best estimate. The uncertainties associated are 0.006 magni- 
tudes given by the spread of our solar analogues divided by the 
square root of number of our solar analogues (namely 10); a 
0.007 magnitudes given by a generous uncertainty in the clus- 
ter metallicity (±0.03). Zero points uncertainties in photometry 
are the dominant source and they account for 0.008. All these 
errors are summed quadratically. To this, the uncertainty in the 
clust er reddening determination, w hich is assumed to be 0.004 
mag jTavlorll2007tlAn et alj|2007l) . is linearly added. 

An additional hidden source of systematic effects, which 
might add a bias towards redder colours, might be present, and 
this is the presence of unidentified binaries. A typical red, faint 
companion will make the stars to appear slightly brighter and 
slightly redder than what they should be, still influencing very 
little the spectroscopic T e ff determination. Given our radial ve- 
locity selection, only a few binaries should be left in our sam- 
ple, and with low mass objects. We cannot quantitatively ac- 
count for their presence, but we shall keep this possibility in 
mind. 

The B - V value found is somewhat in the middle be- 
tw een the majority of the 'old' determinations (see Table 2 



of Barry et al.lll978l which found < (B - V) >=0.667 av- 



eraging most of previous measurements), and the most re- 
cent determinat i ons, w hich, as summarized, for instance, by 
Holmberg et al. (2006), tend to find (B — V) in the range be- 
tween 0.62 and 0.64. None of these results are formally in dis- 
agreement with ours, but we can exclude the estimates at the 
edges of the distribution. 

We think that this estimate is very robust, because our re- 
sults are based on a very few steps and assumptions. We as- 
sume that the metallicity of M67 is essentially solar, and this 
fact is agreed on by all latest works. We determine the T e g in a 
differential way from the Sun, on spectra taken with the same 
instrument, and using two sensitive methods (line-depth ratios 
and Ha wings). We prove that the stars are indeed very close to 
the Sun showing how their spectra overlap with the solar one. 
The stars observed are still on the main sequence. 

6. Cluster distance 

The average V magnitude of the 10 solar analogues 
is 14.583 mag, which must be corrected for reddening: 
Ay=3. 1x0. 041=0. 127, implying a de-reddened V magnitude 
of 14.456. With a solar absolute V magnitude of 4.81 



parison between his colours and those of Montgomery et all dBessell et al.lll998l) . the distance modulus of M67 is of 9.65. 



( 1993), finding an overall zero point shift in B - V of 8 mmag 
(the Sandquist's B - V are bluer than the Montgomery ones). 
The fact that the agreement for these 10 solar stars is better 
than this systematic shift might be due to a statistical fluctu- 
ation, it is on the other hand quite common that calibrations 
agree at best for solar stars. We will nevertheless consider a 
0.008 uncertainty in the colour as introduced by the adopted 
photometry. 

The simple average of the most recent estimate of the M67 
metallicity gives [Fe/H]=0.01. This would require a correction 
of 2 mmag towards the blue for the solar colors derived from 



As mentioned in the previous Section, a correction might be 
needed if the metallicity differs substantially from the solar 
one (of up to 0.08 magnitudes for [Fe/H]=0.05, but we con- 
sider such a large difference in metallicity very unlikely). A 
correction of 0.002 mag, corresponding to a metallicity of 
[Fe/H]=0.01 (used in the previous Section) would bring to a 
distance modulus of 9.63. 

This determinati on is in excellen t agreement with two re- 
cent determinations: I An et alJ (120071) who e stimate a distance 
modulus of 9.61 and ISandquist et alj d2004 . who find 9.60, 
both using the same reddening we adopted. 



L. Pasquini et al.: True solar analogues in M67 



9 



The associated error given by the spread around the aver- 
age magnitude is of 0.060 (i.e. 0.19/ VlO) magnitudes. Other 
sources of uncertainty in the distance modulus will be given by 
the error in reddening, which accounts for 0.012 magnitudes, 
and by the uncertainty on [Fe/H]. If we assume an error on 
[Fe/H] of 0.03 dex, this accounts for 0.05 magnitudes in the 
distance modulus. Summarizing, our best estimate of the dis- 
tance modulus is: 9.63±0.06 sta t ± 0.05 sys . 

A full comparison of our estimate with those present in 
literature is beyond the scope of this work. We find however 
remarkable the agreement betwee n ou r distance esti mate and 
the ones of ISandquist et al.1 (2004) and I An et al.1 ( 12007b in par- 
ticular when considering that our method is independent with 
respect to theirs. 

7. Conclusion 

By using selected observations with FLAMES/GIRAFFE at the 
VLT, we have made a convincing case that the open cluster 
M67 hosts a number of interesting potential solar twins, and 
we have identified them. We have computed spectroscopic ac- 
curate effective temperatures for all the stars with two methods. 
The color-temperature relationships we derive can be used to 
determine temperatures for MS solar- me tallicity stars. 

By computing the average solar twin colours, we have ob- 
tained a precise estimate of the solar (B - V): (B - V) = 
0.649 + 0.016. 

By averaging the magnitude of the solar twins, we have de- 
termined an accurate distance modulus for M67: 9.63+0.06 st at± 
0.05 sys , which is in excellent agreement with the most recent 
estimates, which were based on different, independent meth- 
ods and data sets. 

We have determined for all the stars Li abundances, con- 
firming the presence of a large Li spread among the solar stars 
of this cluster, but showing for the first time, that the Li extra- 
depletion appears only in stars cooler than 6000 K. The candi- 
date solar twins have Li abundance similar to that of our star, 
indicating that they also share with the Sun a similar mixing 
history. 
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Table A.l. Objec t ID, coordinates (Equinox J2000, Epoch J2000.13), photometry, and proper mo tions of the targets (see 



Yadav et al.l 120081 for details). In the last column the name of the stars according to Sanders] d 19771) or to Montgomery et al 
( 1993b aregiven. For a few stars two names are given, because both stars are within 1 arc sec circle, according to SIMBAD. 



Object 


a 




5 


B±AB 


V±AV 


I ±AI 


M* ± A Mx 


My * A/ij 


Name 




O 




O 


(mag) 


(mag) 


(mag) 


(mas/yr) 


(mas/yr) 




Obj219 


132.893537 


1 1 


632623 


13.715±0.009 


13.100±0.008 


12.425±0.005 


1.96*1.31 


8.63±3.03 


SI 197 


Obj266 


132.860339 


1 1 


643584 


14.212±0.006 


13.601±0.006 


12.939±0.007 


-0.24±2.38 


0.18±3.45 


S944 


0bj285 


132.849450 


11 


647816 


15.165±0.006 


14.461±0.00() 


13.713±0.002 


-0.42±1.78 


-2.80±1.37 


S945 


0bj288 


132.900657 


11 


648980 


14.494±0.013 


13.857±0.004 


13.160±0.005 


0.00±0.77 


1.61 ±2.44 


S1201 


Obj291 


132.937203 


1 1 


649712 


14.090±0.010 


13.478±0.009 


12.807±0.018 


0.18±1.13 


7.26±3.03 


SI 202 


Obj342 


132.823311 


1 1 


660013 


14.118±0.004 


13.503±0.005 


12.841±0.008 


0.77±2.14 


5.59±4.22 


S948 


0bj349 


132.979809 


1 1 


661163 


14.978±0.006 


14.301±0.002 


13.614±0.007 


2.92±8.87 


2.92±4.82 


MMJ6241/S1423 


Obj350 


132.836621 


1 1 


661145 


14.226±0.005 


13.624±0.010 


12.955±0.006 


-1.96±1.07 


0.18±1.49 


S950 


0bj364 


132.794912 


11 


664138 


15.288±0.005 


14.584±0.010 


13.809±0.006 


-1.37±1.55 


0.24±1.43 


S951 


Obj401 


132.829348 


1 1 


671034 


14.268±0.006 


13.661±0.007 


13.009±0.003 


-3.81±2.56 


1.31±1.96 


S954 


Obj437 


132.827827 


1 1 


676879 


14.600±0.006 


13.998±0.008 


13.331±0.002 


2.98±0.71 


-2.02±3.69 


S956 


Obj455 


132.971225 


1 1 


681574 


13.930±0.007 


13.301±0.004 


12.608±0.018 


5.30±5.65 


5.77±2.80 


S1428 


Obj473 


132.809740 


1 1 


685891 


15.142±0.011 


14.443±0.008 


13.731±0.004 


1.31±1.13 


-1.55±3.57 


S958 


Obj571 


132.951947 


1 1 


706357 


14.686±0.012 


14.021±0.014 


13.309±0.004 


-0.89±1.31 


-0.30±3.21 


S1211 


Obj574 


133.024207 


1 1 


706858 


14.309±0.016 


13.653±0.007 


12.932±0.005 


-1.73±1.01 


-0.24±1.07 


MMJ6362/S1431 


Obj587 


132.911415 


1 1 


710387 


14.753±0.006 


14.107±0.006 


13.405±0.020 


-0.89±0.83 


1.31±3.15 


S1213 


Obj613 


132.825417 


11 


715202 


13.907±0.015 


13.254±0.006 


12.594±0.022 


1.78±2.62 


-0.89±3.69 


S964 


Obj637 


132.840991 


11 


721590 


15.191±0.017 


14.489±0.010 


13.751±0.011 


0.65±0.54 


-0.42±2.50 


S966 


Obj673 


132.722266 


11 


727791 


15.062±0.002 


14.356±0.009 


13.569±0.001 


0.00±0.77 


0.24±2.86 


S746 


Obj681 


132.773982 


11 


729705 


14.715±0.006 


14.018±0.016 


13.254±0.008 


1.19±1.31 


1.49±2.38 


S747 


Obj689 


132.875550 


1 1 


730521 


13.783±0.012 


13.120±0.011 


12.436±0.010 


-3.21±2.86 


-0.36±3.87 


S1219 


Obj713 


132.840699 


11 


734734 


14.886±0.017 


14.172±0.009 


13.421±0.017 


1.55±1.37 


-0.48±1.31 


S969 


Obj750 


132.722677 


11 


742964 


14.215±0.007 


13.576±0.014 


12.865±0.005 


0.18±2.08 


1.07±2.20 


S750 


Obj756 


132.947493 


1 1 


745099 


15.393±0.011 


14.695±0.008 


13.923±0.009 


3.03±0.89 


5.06±1.73 


S1222 


Obj769 


132.959247 


11 


749185 


14.119±0.004 


13.478±0.003 


12.771±0.005 


-0.36±1.13 


0.12±2.98 


S 1224a 


Obj778 


132.836684 


1 1 


750684 


13.716±0.010 


13.093±0.011 


12.41 1±0.011 


1.37±1.78 


2.98±4.40 


S976 


Obj809 


132.758539 


11 


755306 


15.696±0.007 


14.959±0.015 


14.162±0.004 


-1.13±1.61 


-1.49±1.84 


S754 


Obj851 


132.854122 


11 


761931 


14.730±0.007 


14.113±0.004 


13.449±0.007 


0.71 ±3.09 


2.98±1.07 


S982 


Obj88() 


132.770112 


11 


765793 


14.202±0.009 


13.544±0.014 


12.841±0.002 


-0.30±1.07 


-1.07±1.61 


S757 


Obj905 


132.746798 


1 1 


770272 


14.038±0.004 


13.434±0.010 


12.748±0.015 


1.01±1.73 


-1.73±2.74 


S758 


Obj911 


132.791312 


1 1 


771367 


15.220±0.006 


14.547±0.010 


13.785±0.010 


0.36±1.31 


1.19±1.31 


S991 


Obj917 


133.021288 


1 1 


772611 


15.497±0.007 


14.755±0.008 


13.923±0.006 


-0.54*2.56 


2.44±2.62 


S1442 


Obj971 


132.884919 


1 1 


779352 


15.287±0.005 


14.592±0.004 


13.793±0.000 


-0.12±0.59 


O.36±0.3O 


S1246 


Obj986 


132.905212 


1 1 


782141 


14.646±0.0()8 


14.007±0.003 


13.283±0.001 


0.24±1.25 


1.19±0.54 


S1247 


Obj988 


132.892171 


11 


782156 


14.819±0.005 


14.180±0.001 


13.475±0.002 


-0.65±0.65 


-0.24±1.49 


S1248 


ObjlOlO 


132.850486 


1 1 


785927 


14.104±0.007 


13.478±0.009 


12.781±0.021 


-0.36±0.95 


3.57±1.43 


S2209 


Objl()32 


132.985778 


11 


790265 


14.997±0.005 


14.358±0.003 


13.649±0.001 


3.33±3.33 


-0.24±0.48 


S1449 


Objl036 


132.851246 


1 1 


791211 


15.678±0.018 


14.947±0.003 


14.164±0.005 


-0.36±0.24 


-1.01±1.25 


S1004 


Objl 051 


132.922899 


1 1 


793371 


14.726±0.007 


14.09()±0.004 


13.382±0.004 


-1.55±0.54 


1.96±1.07 


S1252 


Objl062 


132.900316 


1 1 


796392 


15.144±0.001 


14.477±0.008 


13.745±0.004 


1.25±1.07 


-0.12±2.68 


S1255 


Objl()67 


133.014592 


1 1 


796693 


15.201±0.005 


14.559±0.009 


13.824±0.002 


0.24±2.08 


-0.30±1.25 


S1452 


Objl()75 


132.912716 


1 1 


798715 


14.386±0.000 


13.712±0.006 


12.992±0.004 


1.19±1.96 


-2.80±3.33 


S1256 


Objl()88 


132.953942 


1 1 


800602 


15.151±0.004 


14.492±0.001 


13.760±0.007 


0.77±1.90 


-3.99±2.56 


S1258 


Objl090 


132.869575 


11 


800607 


14.450±0.004 


13.800±0.005 


13.040±0.010 


0.18±1.73 


-1.55±1.07 


S1011 


Objl091 


132.853518 


11 


801110 


15.237±0.011 


14.513±0.007 


13.669±0.008 


0.12±0.83 


-0.48±1.55 


S1012 


ObjllOl 


133.035214 


1 1 


802908 


15.377±0.001 


14.675±0.004 


13.903±0.001 


-0.06±1.31 


-2.08±2.50 


MMJ6384 


Objll()8 


132.855293 


11 


803762 


14.878±0.005 


14.177±0.001 


13.386±0.006 


-0.30±0.95 


-0.42±0.71 


S1014 


Objl 129 


132.885730 


1 1 


806694 


14.795±0.005 


14.171±0.006 


13.482±0.002 


1.13±0.30 


-1.55*1.31 


S1260 


Objl 137 


132.800104 


1 1 


807413 


15.571±0.002 


14.873±0.008 


14.107±0.009 


-1.90±0.89 


0.83*1.19 


S2213 


Objl 161 


132.981615 


11 


810602 


14.549±0.004 


13.883±0.010 


13.149±0.006 


0.89±4.46 


-6.84±2.86 


S1457 


Objl 163 


132.787098 


1 1 


810483 


14.688±0.009 


13.992±0.004 


13.221±0.002 


0.42±0.65 


8.03*1.19 


S1022 


Objl 194 


132.753356 


1 1 


814655 


15.281±0.002 


14.614±0.010 


13.876±0.002 


0.30±1.01 


-0.42*0.65 


S770 


Objl 197 


132.877372 


11 


815215 


13.921±0.006 


13.315±0.003 


12.618±0.024 


-0.18±1.07 


1.78*0.77 


MMJ5882/S 1264b 


Objl 247 


132.812708 


1 1 


822522 


14.753±0.008 


14.144±0.008 


13.470±0.005 


0.65±1.25 


-2.80*1.96 


S1033 


Objl 303 


132.736097 


11 


831844 


15.318±0.008 


14.641±0.008 


13.899±0.009 


0.06±2.26 


0.89*1.25 


S779 


Objl 304 


132.858515 


1 1 


832075 


15.454±0.015 


14.731 ±0.009 


13.916±0.006 


-0.12±1.31 


-0.12*1.90 


S1041 


0bjl315 


132.994897 


11 


834025 


14.990±0.013 


14.297±0.01 1 


13.544±0.008 


-1.19±4.40 


0.48*2.20 


MMJ6306/S1462 


Objl 3 34 


132.866453 


1 1 


836636 


15.083±0.007 


14.403±0.007 


13.669±0.000 


0.77±1.61 


-1.84*0.89 


SI 048 


Objl 342 


132.826093 


1 1 


838787 


14.935±0.006 


14.285±0.005 


13.547±0.001 


1.37±1.49 


-0.36*0.06 


S1050 


0bjl387 


132.874854 


11 


852505 


14.724±0.004 


14.098±0.002 


13.398±0.000 


0.59±3.75 


-1.55*7.97 


S1283 


Objl 392 


132.749200 


1 1 


853500 


15.527±0.002 


14.81 1±0.004 


14.047±0.000 


-1.49±2.08 


-1.61*0.77 


S785 


Objl 397 


132.882992 


11 


854616 


14.631±0.001 


14.009±0.003 


13.304±0.001 


1.84±2.32 


-3.03*5.89 


S1287 


Objl 424 


132.890226 


11 


862493 


13.825±0.009 


13.203±0.005 


12.501±0.000 


1.49±2.50 


-1.25*3.87 


Check 


0bjl458 


132.762475 


1 1 


873808 


15.716±0.010 


14.977±0.005 


14.186±0.004 


1.01±1.84 


0.12*1.49 


S795 


Objl480 


132.916988 


11 


878716 


14.383±0.008 


13.783±0.008 


13.125±0.002 


-3.99±2.50 


1.25*2.32 


SI 300 


0bjl496 


132.781325 


1 1 


882262 


14.486±0.013 


13.879±0.004 


13.214±0.002 


-0.95±0.12 


-0.30*0.65 


S2224 


Objl 504 


132.864557 


1 1 


884028 


14.796±0.001 


14.171±0.011 


13.474±0.000 


-1.13±0.95 


1.31*1.31 


S1078 


0bjl514 


132.753181 


11 


886527 


15.498±0.003 


14.777±0.004 


14.008±0.000 


-0.12±1.13 


1.73*1.37 


S802 


0bjl587 


132.846438 


1 1 


901394 


14.804±0.015 


14.163±0.004 


13.469±0.006 


-0.95±1.73 


-1.01*1.78 


S1087 


Objl 622 


132.801216 


1 1 


906389 


14.788±0.004 


14.156±0.002 


13.459±0.004 


1.43±1.25 


0.95*0.65 


S1089 


Objl 680 


132.883591 


1 1 


919069 


14.291±0.015 


13.646±0.006 


12.951±0.004 


1.01*1.13 


0.18*1.07 


S1314 


Objl 706 


133.010317 


1 1 


926166 


15.468±0.005 


14.745 ±0.0 17 


13.981±0.006 


0.54±0.95 


0.30*0.54 


MMJ6341/S1481 


0bjl716 


132.866205 


1 1 


928044 


13.918±0.010 


13.299±0.009 


12.625±0.005 


-1.96±1.84 


4.82*3.15 


S1092 


Objl 722 


132.828014 


1 1 


930478 


14.731 ±0.002 


14.130±0.006 


13.449±0.009 


-3.57±26.24 


-2.44*2.20 


SI 093 


Objl 735 


133.001728 


11 


935288 


14.993±0.012 


14.332±0.010 


13.617±0.007 


1.96±1.78 


-2.92*1.55 


S1483 


Objl 75 8 


132.746824 


11 


943570 


13.860±0.056 


13.207±0.015 


12.545±0.007 


-0.89±0.48 


3.45*1.67 


MMJ5342/S816 


Objl 768 


132.906577 


1 1 


945718 


15.060±0.003 


14.404±0.004 


13.684±0.004 


0.77±1.13 


-1.01*0.89 


MMJ6028/S1320 


Objl 778 


132.737778 


1 1 


947424 


15.679±0.004 


14.948±0.004 


14.155±0.003 


0.83±1.96 


-0.65*0.65 


MMJ5310/S820 
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Table A.2. Radial velocities, effective temperatures, and lithium abundances of the 59 stars retained as possible single members. 
The values in Table are derived by adding the difference of r ei f between the stars and the Sun to the solar temperature (5777 K). 
The solar temperature derived from the GIRAFFE spectrum and our calibrations is of 5792 and 5717 K for the LDR and Ha 
methods, respectively. The 10 best solar twins candidates are indicated in bold face. The S /N ratio/pixel of the co-added spectra 
varies between 80 and 110, depending on the magnitude of the stars. 



Object 


Vrad ± AV„, d 


T Lbk , AT Lbk 

eff efF 


oir 


EW Ll+ F C 


log£W L i 


log/V(Li) L,E 


logiV(Li) ,NL1E 




(kms- 1 ) 


(K) 


(K) 


(mA) 


(mA) 






Sun 




5777±27 


sm 


12.1 


0.3 


0.8 


0.8 


219 


33.85±0.28 


6243±54 


6110 


22.7 


1.2 


2.1 


2.1 


266 


32.95±0.46 


6147±63 


6060 


53.2 


1.7 


2.6 


2.5 


285 


33.72±0.67 


5836±67 


sm 


1.5 


<0.0 


0.6 


<0.6 


288 


34.26±0.40 


6004±67 


6010 


42.4 


1.5 


2.3 


2.3 


291 


32.13±0.29 


6177±57 


6160 


59.6 


1.7 


2.6 


2.6 


349 


34.28±0.42 


5952±78 


5917 


10.4 


<0.0 


0.7 


<0.7 


350 


32.56±0.31 


6024±52 


6010 


68.2 


1.8 


2.6 


2.6 


401 


32.64±0.37 


6165±64 


6110 


44.9 


1.6 


2.5 


2.4 


473 


34.74±0.38 


5919±76 


5807 


2.7 


<0.0 


0.7 


<0.7 


587 


33.06±0.39 


6077±65 


6060 


37.2 


1.4 


2.3 


2.3 


613 


33.05±0.29 


6202±45 


6110 


61.3 


1.7 


2.7 


2.6 


637 


34.00±0.50 


5806±65 


sin 


17.8 


0.9 


1.4 


1.4 


673 


32.44±0.57 


5639±63 


5747 


19.1 


0.9 


1.4 


1.4 


689 


32.88±0.24 


6093±41 


6110 


42.3 


1.5 


2.4 


2.3 


750 


33.13±0.28 


5918±48 


5927 


16.4 


0.9 


1.5 


1.5 


769 


34.47±0.29 


5984±46 


6010 


29.4 


1.3 


2.1 


2.0 


778 


33.45±0.24 


6114±39 


6060 


15.6 


0.8 


1.7 


1.7 


809 


31.87±0.48 


5667±78 


5537 


6.9 


<0.0 


0.4 


<0.5 


851 


33.47±0.31 


5948±61 


6060 


12.9 


0.6 


1.3 


1.3 


911 


32.08±0.35 


5885±67 


5837 


13.4 


0.6 


1.2 


1.2 


988 


32.07±0.28 


5935±53 


6060 


28.2 


1.3 


2.0 


2.0 


1032 


34.02±0.47 


5955±60 


6010 


17.8 


0.9 


1.6 


1.6 


1036 


33.39±0.48 


5612±66 


5537 


20.9 


1.0 


1.4 


1.4 


1051 


32.21±0.29 


6081 ±55 


6060 


34.7 


1.4 


2.2 


2.2 


1062 


32.64±0.45 


5926±55 


5867 


21.0 


1.1 


1.7 


1.7 


1067 


33.37±0.35 


5929±69 


5917 


11.4 


0.4 


1.0 


1.0 


1075 


33.13±0.24 


5871 ±48 


5917 


10.7 


0.0 


0.6 


0.7 


1088 


32.87±0.29 


5890±59 


5867 


9.3 


<0.0 


0.6 


<0.7 


1090 


33.34±0.32 


6086±54 


6010 


41.2 


1.5 


2.3 


2.3 


1101 


32.72±0.34 


5756±60 


5717 


6.9 


<0.0 


0.5 


<0.6 


1129 


33.92±0.30 


5959±51 


6010 


32.4 


1.4 


2.1 


2.1 


1137 


33.59±0.48 


5741 ±69 


5627 


5.3 


<0.0 


0.5 


<0.6 


1194 


33.30±0.40 


5766±64 


5837 


5.3 


<0.0 


0.5 


<0.6 


1197 


34.26±0.28 


6207±44 


6110 


28.1 


1.3 


2.2 


2.2 


1247 


32.44±0.51 


5994±60 


6010 


30.2 


1.3 


2.1 


2.1 


1303 


32.65±0.41 


5716±64 


5717 


10.1 


<0.0 


0.5 


<0.6 


1304 


33.60±0.39 


5704±64 


5717 


7.9 


<0.0 


0.5 


<0.5 


1315 


32.55±0.34 


5874±58 


5867 


15.6 


0.8 


1.4 


1.4 


1334 


32.37±0.45 


5957±57 


5957 


29.2 


1.3 


2.0 


2.0 


1387 


33.35±0.24 


6090±58 


6060 


37.5 


1.5 


2.3 


2.3 


1392 


33.68±0.57 


5716±63 


5687 


6.3 


<0.0 


0.5 


<0.6 


1458 


32.55±0.56 


5640±65 


5567 


6.1 


<0.0 


0.4 


<0.5 


1496 


34.22±0.48 


6173±54 


6160 


63.1 


1.7 


2.7 


2.6 


1504 


33.39±0.48 


5934±56 


6060 


40.7 


1.5 


2.2 


2.2 


1514 


33.33±0.49 


5613±67 


5597 


25.6 


1.2 


1.6 


1.6 


1587 


32.38±0.62 


5975±55 


6010 


28.6 


1.3 


2.0 


2.0 


1622 


33.26±0.61 


6043±57 


6010 


34.6 


1.4 


2.2 


2.2 


1716 


33.84±0.59 


6030±40 


6060 


40.8 


1.5 


2.3 


2.3 


1722 


33.87±0.52 


6007±62 


6010 


38.1 


1.5 


2.3 


2.2 


1735 


33.18±0.59 


5959±59 


5960 


10.9 


0.2 


0.9 


0.9 


1758 


33.78±0.67 


6221 ±52 


6160 


36.3 


1.4 


2.4 


2.4 


1768 


33.98±0.57 


5844±57 


5927 


40.1 


1.5 


2.1 


2.1 


1787 


33.44±0.57 


5768±70 


5807 


12.9 


0.5 


1.1 


1.0 


1788 


33.38±0.80 


5886±60 


5867 


23.4 


1.1 


1.8 


1.8 


1852 


32.30±0.40 


6009±63 


6010 


32.2 


1.4 


2.1 


2.1 


1903 


32.76±0.42 


5609±72 


5687 


6.8 


<0.0 


0.4 


<0.5 


1948 


32.86±0.29 


6164±63 


6010 


43.0 


1.5 


2.4 


2.4 


1955 


32.63±0.60 


5961 ±76 


5837 


35.8 


1.4 


2.2 


2.2 


2018 


31.78±0.43 


5693±74 


5777 


8.7 


<0.0 


0.4 


<0.5 
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Table A.3. Radial velocities of likely binaries or non members. For these stars each of the three RV measurements is given. Fields 
with blanks indicate that problems were present with the cross-correlation profile of these objects. 



Object 




AV 1 , 

rad 






vL 




342 


17.74 


0.72 


57.20 


1.45 


40.75 


0.73 


364 


36.93 


0.43 


36.58 


1.18 


36.48 


0.74 


437 


45.87 


0.65 


45.54 


0.83 


45.68 


0.52 


455 


38.07 


0.50 


14.62 


1.24 


33.41 


0.47 


571 


31.73 


0.49 


31.07 


0.66 


31.47 


0.56 


574 


42.88 


0.46 


42.91 


0.85 


42.08 


0.62 


681 


24.14 


0.80 


56.70 


1.22 


59.70 


0.85 


713 


35.34 


0.46 


35.03 


0.64 


34.91 


0.33 


756 


46.24 


0.50 


45.54 


0.77 


45.47 


0.41 


880 


39.62 


6.67 


/ 


/ 


33.55 


5.81 


905 


34.78 


0.35 


37.30 


1.02 


74.65 


0.89 


917 


18.66 


0.86 


15.59 


1.36 


/ 


/ 


971 


36.13 


0.51 


36.13 


0.70 


36.22 


0.63 


986 


16.54 


0.53 


27.69 


0.84 


49.08 


0.68 


1010 


17.22 


0.97 


30.77 


0.63 


34.29 


0.75 


1091 


27.59 


3.27 


27.03 


0.62 


25.77 


0.78 


1108 


15.92 


1.00 


/ 


/ 


/ 


/ 


1161 


35.53 


0.65 


35.49 


1.21 


22.45 


1.55 


1163 


10.10 


0.66 


65.99 


1.32 


70.80 


1.24 


1342 


37.97 


0.69 


28.40 


1.06 


34.93 


0.65 


1397 


35.43 


1.12 


34.91 


0.85 


34.88 


0.54 


1424 


54.44 


1.62 


57.22 


1.31 


56.04 


1.58 


1480 


22.86 


1.02 


25.56 


1.09 


28.27 


0.48 


1680 


37.65 


1.29 


35.65 


0.96 


32.27 


0.81 


1706 


35.79 


1.05 


35.40 


0.90 


35.77 


0.61 


1778 


29.59 


1.30 


29.92 


2.96 


31.19 


1.26 


1842 


30.09 


0.69 


29.88 


1.20 


30.52 


0.47 


1862 


31.31 


0.63 


30.29 


2.03 


31.27 


0.43 


1957 


21.17 


0.52 


20.15 


0.83 


21.21 


0.57 


2016 


21.64 


0.46 


25.96 


0.94 


65.63 


0.66 


2017 


126.58 


0.84 


125.09 


1.62 


126.01 


0.63 



